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ABSTRACT 

We study properties of luminous X-ray pulsars using a simplified model of the accretion 
column. The maximal possible luminosity is calculated as a function of the neutron 
star (NS) magnetic field and spin period. It is shown that the luminosity can reach 
values of the order of 10 4o ergs -1 for the magnetar-likc magnetic field (B )> 10 14 G) 
and long spin periods (P )> 1.5 s). The relative narrowness of an area of feasible 
NS parameters which are able to provide higher luminosities leads to the conclusion 
that L ~ 10 4o ergs _1 is a good estimate for the limiting accretion luminosity of a 
NS. Because this luminosity coincides with the cut-off observed in the high mass X- 
ray binaries luminosity function which otherwise does not show any features at lower 
luminosities, we can conclude that a substantial part of ultra-luminous X-ray sources 
are accreting neutron stars in binary systems. 
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1 INTRODUCTION 

Accretion of matter on to a compact object in a binary sys¬ 
tems is an extremely efficient way of producing high-energy 
photons. However, the luminosity that can be reached in this 
process is limited by different factors. In the case of spheri¬ 
cal accretion, the luminosity is restricted by the well known 
Eddington luminosity Z/Edd = 4,nGMmnc/[aT{l + X)] ~ 

2 x 10 38 ergs -1 , for a typical mass of the neutron star (NS) 
M = 1.4Mq (here rriH is hydrogen mass and X is its mass 
fraction). On the other hand, the luminosity of highly mag¬ 
netized NS can exceed Z/Edd substantially for two reasons: 
the interaction cross-section can be much below the Thom¬ 
son value and the radiation can escape in the direction per¬ 
pendicular to the incoming flow of the material which is 
funnelled to the accretion column (Basko & Sunyaev 1976; 
Lyubarskii & Syunyaev 1988). The radiation pressure is bal¬ 
anced here by strong magnetic pressure instead of gravity 
and the corresponding flux can exceed the Eddington flux 
significantly. 

Some transient X-ray pulsars indeed produce 
super-Eddington luminosities (Lucke et al. 1976; 
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Kouveliotou et al. 1996; Tsygankov et al. 2006). Fur¬ 
thermore, a recent revolutionary discovery by the NuSTAR 
observatory of pulsations from an ultra-luminous X-ray 
source (ULX) X-2 in galaxy M82 (Bachetti et al. 2014) 
implies that accreting NS can reach luminosities of about 
10 4 °ergs -1 , two orders of magnitude above the Eddington 
limit. Until recently one thought that either the intermediate 
mass black holes 10 2 — 10 5 Mq (Colbert & Mushotzky 1999) 
or the super-critically accreting stellar mass black holes (e.g. 
Poutanen et al. 2007) are the most probable sources of such 
high X-ray luminosities (see review by Feng & Soria 2011). 
Although highly-magnetized NSs have previously been 
suggested to power ULXs (Medvedev & Poutanen 2013), it 
was rotation rather than accretion that was responsible for 
high energetics. 

From the standard theory of accretion onto magne¬ 
tized NSs it is known that the configuration of emit¬ 
ting region highly depends on the mass accretion rate 
(Basko & Sunyaev 1976). Particularly, at low luminosity the 
X-ray emission originates from the hot spots at the NS sur¬ 
face. As soon as the pulsar luminosity exceeds the criti¬ 
cal value (~ 10 37 ergs -1 ) an accretion column begins to 
rise above the NS (Basko & Sunyaev 1976; Mushtukov et al. 
2015) and provides a possibility for the luminosity to exceed 


©2015 The Authors 


2 A. A. Mushtukov et al. 


the Eddington limit. The radiation pressure is balanced by 
the magnetic pressure in the accretion column, which sup¬ 
ports the column. The height of the accretion column in¬ 
creases with the mass accretion rate and can be as high as 
the NS radius (Lyubarskii & Syunyaev 1988). Even higher 
accretion columns are probably possible but their luminosity 
are not expected to grow significantly because of a rapidly 
decreasing accretion efficiency. 

In this paper, we develop a simplified model of ac¬ 
cretion onto magnetized NSs based on the suggestions 
by Basko & Sunyaev (1976) and Lyubarskii & Syunyaev 
(1988). For the first time, we take into account exact Comp¬ 
ton scattering cross-sections in strong magnetic field. The 
structure and the maximal luminosity of the accretion col¬ 
umn are calculated. We determine an upper limit on the NS 
luminosity to be close to 10 4O ergs -1 , which coincides with 
the position of a sharp cut-off in the luminosity function 
of high-mass X-ray binaries (Mineo, Gilfanov, & Sunyaev 
2012) indicating that a large fraction of bright X-ray bi¬ 
naries are likely accreting magnetized NSs. A possibility of 
extremely high mass accretion rate (M J> 10 19 gs -1 ) onto a 
magnetized NS is discussed. We find that such high M imply 
extremely high, magnetar-like magnetic field strength. 

2 ACCRETION COLUMN MODEL 
2.1 Analytical estimates 

Let us consider an X-ray pulsar accreting at a suffi¬ 
ciently high rate. If luminosity exceeds the critical value 
L* (Basko & Sunyaev 1976; Mushtukov et al. 2015), the 
radiation-dominated shock is formed in the accreting gas. 
After the shock, the hot gas settles to the NS surface in the 
magnetically confined accretion column, whose height de¬ 
pends on the accretion rate and the magnetic field strength 
(see Fig. 1). At high accretion rates that are of interest here, 
the optical depth of the settling flow along and across the 
magnetic field direction is significantly larger than unity. 

Our model is based on a description of optically 
thick accretion column proposed by Lyubarskii & Syunyaev 
(1988). The cross-section of the accretion column is a thin 
annular arc of width do which is much less than its length 
lo. Then the area of the accretion column base is Sd = lodo- 

The gravitational force, which acts on the gas in the 
accretion column, will be offset by the radiation pressure 
gradient only. 1 The gas pressure becomes important at 
virial temperature ~ 10 12 K at the NS surface, which 
cannot be reached due to effective electron-positron pair 
creation and further neutrino cooling at much lower tem¬ 
peratures ~ 10 10 K (Beaudet, Petrosian & Salpeter 1967; 
Basko & Sunyaev 1976). The radiation pressure gradient at 
the centre of the accretion channel can be presented as fol¬ 
lows: 

dP Ta d(x = 0, h) GM 

dh P (R + h) 2 ’ U 

1 Photon diffusion time from the centre of the column is tdiff — 
dor/(2c) ~ 10“ 5 s. The matter settling time in the accretion col¬ 
umn is tdyn > 7H/vff(H) ~7x 10“ 4 s (for the accretion column 
of height H = R). As a result, tdiff tdyn and one can neglect 
in the first approximation the effects of plasma movement. 



Figure 1. The structure of the accretion column. At the NS 
surface the cross section of the accretion channel is a thin annular 
arc of length Iq and thickness do. The position of the radiation- 
dominated shock varies inside the accretion channel, because the 
radiation energy density drops sharply towards the column edges. 
As a result the accretion column (i.e. the settling flow) height 
reaches its maximum inside the channel and decreases towards 
its borders. 

where M and R are NS mass and radius, h is the height 
above the stellar surface, p is the local mass density in the 
column. This equation gives an upper limit of P r a d, since we 
ignore the contribution of the gas pressure and deviations of 
the hydrostatic equilibrium due to a subsonic settlement of 
the gas in the column. This approximation is equivalent to 
a suggestion that the radial radiation flux at every height of 
the accretion column equals the local Eddington flux 

E Edd = — GA ( (2) 

kII [R + h) 2 

where «|| is the opacity along the magnetic field lines. As¬ 
suming constant density and taking a zero boundary con¬ 
dition ( P T s.d(x = 0 , H) = 0) for the hydrostatic equilibrium 
equation, we can get a simple analytical expression for the 
radiation pressure distribution along the column height 

„ , , N GM H-h 

P^d{x - 0, h) » (3) 

where H is the height of the shock. Assuming further the 
optically thick column, we can estimate the flux emergent 
from the column side at a given height h using a plane- 
parallel diffusion approximation for the radiation transfer 
across the accretion channel: 

dPrad^x, h') Fj_(x 1 h') Fj_ esc (h) 2x , 

ax c c do 

where a: is a distance from the centre of the accretion channel 
(see Fig.l), k± is the opacity across the B-field lines, and 
F±,esc(h) = F±(do/2, h) is the escaping flux. The second 
equality here was obtained from the condition Fj_ (0, h) = 0 
and assuming that the flux increases linearly towards the 
edges of the column (where x = do/2). Then the solution of 


MNRAS 000, 1-11 (2015) 






















The maximum luminosity for magnetized neutron stars 3 


equation (4) is 


2.1.1 Consistency check 


-Trad(o:, h) 


r 0 F± y e SC (h) ( 4a; 2 \ 

4c \f dl ) 


(5) 


where to = p K,±do is the optical depth across the accretion 
channel. This solution allows us to find an emergent flux: 


Fj_,eac (h) 


4c P ra d (x = 0, h) 

TO 


4c GM H-h 

ki d 0 R + H R + h’ ( ' 


which is much higher than the Eddington flux for the case 
of tall (H do) accretion columns 


Tj_ J esc(h) 



H-h R+h 
do R + H 


^Edd- 


(7) 


The total accretion luminosity of a NS can be estimated by 
integrating the emergent flux of the surface of two columns: 


L = 4Z 0 [ Tx,esc(h) dh 

Jo 


16 ^ 



^GM 

K± 


« 38 


/ Ip/dp 

V 50 



( 8 ) 


where kt ~ 0.34 cm 2 g 1 is the Thomson electron scattering 
opacity of the solar mix plasma and 


/ 



■4+f) 


H 

R + H ' 


(9) 


The function f(x) grows quadratically f(H/R) ~ ( H/R ) 2 
at small column heights H R and grows logarithmically 
at H > R. For the simple analytical estimate (8) we as¬ 
sumed that the optical thickness across the column is con¬ 
stant. The dependence of the optical thickness on the height 
will be taken into account accurately in our numerical model 
(see Section 2.2). The maximum X-ray pulsar luminosity 
can be estimated from equation (8) substituting H = R 
(Basko & Sunyaev 1976) 

( 10 ) 

We have suggested earlier that the accretion column 
thickness does not depend on the height above the surface. 
In fact, the radiation pressure drops down towards the ac¬ 
cretion channel edges (see equation (5)). As a result, the ra¬ 
diation pressure cannot support the same column height all 
over the accretion channel. The height has its maximum in 
the middle of the channel and decreases towards its edges. 
Therefore, the column thickness has to decrease with the 
height as well (see Fig. 1). Taking the radiation pressure 
P TS .d(x,h = 0) at the bottom of the accretion column we 
can evaluate the column height H x (i.e. the position of the 
radiative shock) as for a given x using equation (3) where H 
is replaced with H x . This height drops from the maximum 
value H at the centre to the zero value at the column edge: 

H x H f 4a; 2 

R + H x R + H \ d\ 

Therefore, the radial cross section of homogeneous column 
of a small height (H R) has a roughly parabolic form 

(Fig. 1). 



Let us estimate the values of the basic parameters in or¬ 
der to justify the proposed model. Let us consider a NS 
of mass M = 1.4Mq and radius R = 10 km accreting at 
a rate M « 5 x 10 18 gs -1 (which gives accretion lumi¬ 
nosity L = 10 39 ergs -1 ). We take the column base area 
5 d = 5 x 10 9 cm 2 and its geometrical size are lo = 5 x 10 5 cm 
and do = 10 4 cm, which are expected for a given mass 
accretion rate for standard magnetic field of B = 10 12 G 
(Basko & Sunyaev 1976). We also assume here k± = kt- 

According to equation (8) the corresponding height of 
the accretion column is H ~ 0.75 R. The free-fall velocity at 
the height H above the NS surface vg = \J2GM / (R + H) m 
1.5 x 10 10 cms -1 and the mass density for given mass ac¬ 
cretion rate pg ~ 3.5 x 10 -2 gem” 3 . Because the mat¬ 
ter is decelerated by the radiation-dominated shock at 
the top of the accretion column (Basko & Sunyaev 1976), 
the plasma velocity falls down to v = vg/7 (see, e.g. 
Lyubarskii & Syunyaev 1982; Becker 1998), and the mass 
density goes up to p ~ 0.245 gem -3 . As a result the 
Thomson optical thickness across the accretion channel is 
ro = pKxdo ~ 500, which means that the accretion column 
is optically thick as it was assumed earlier. 

The radiation pressure and the emergent side flux given 
by equations (3) and (6), respectively, near the NS sur¬ 
face are P ra d ~ 3.5 x 10 19 dyncm -2 and F± a 8 x 
10 28 ergs -1 cm -2 . Because the column is optically thick we 
can assume thermodynamical equilibrium deep inside it and 
evaluate the plasma temperature from the radiation field 
energy density: 


Trad 


£rad 

“3” 


aT 4 
3 ’ 


( 12 ) 


where a is the radiation constant. This relation gives T « 
3.4 x 10 s K or kT pc 30 keV. Therefore, the gas pressure 
P g = pkT/(p s mn) ~ 5 x 10 15 dyncm -2 is much smaller 
than the radiation pressure -Prad, as we assumed earlier (here 
p,s = 0.62 is the mean particle weight for fully ionized solar 
mix plasma). The effective temperature corresponding to the 
escaping flux is determined by the relation F± = asBlj 1 is 
Tr s; 1.9 x 10 8 K or kT± « 16 keV. This value is close 
to the cut-off energy observed in the spectra of bright X- 
ray pulsars, especially if the gravitational redshift is taken 
into account (Filippova et al. 2005). We see thus that the 
assumption made constructing the accretion column model 
are reasonable and the model is self-consistent. We consider 
this analytical model as a base for further numerical studies. 

In the analytical estimates, we have considered a homo¬ 
geneous column with a constant cross-section area. For the 
numerical model described in the next section, we account 
for the column inhomogeneity and take the exact opacity of 
highly-magnetized plasma. 


2.2 Geometry and structure of the accretion 
column 

Using the approach described in Section 2.1, we can con¬ 
struct a numerical model of the accretion column. Global 
input model parameters are the NS mass M and radius R, 
the mass accretion rate M, and the polar magnetic field 
strength B. The magnetic field is assumed to be dipole with 
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the magnetic field strength decreasing with distance as 

Bh ~ B {^ir) =BS ^’ (13) 

where Sh is the cross-section area of the accretion channel at 
a given height h above the stellar surface. The width and the 
length of the accretion channel depends on height as follows: 

-m" <«> 


The geometrical parameters of accretion column base lo and 
do can be derived from the global parameters of the model 
(see Sect. 2.4). It has to be mentioned that the thickness of 
a sinking region of the accretion column dh at some non¬ 
zero height h is always smaller than its base thickness do 
because of the parabolic shape of the shock surface. The 
column thus is surrounded by almost free-fall accreting mat¬ 
ter of geometrical thickness ~ (d — dh)/ 2, which affects the 
distribution of the emergent radiation over the directions 
(Lyubarskii & Syunyaev 1988; Poutanen et al. 2013). 

We consider a two-dimensional radial cross section of 
the column using two geometrical coordinates h and x 
(Fig. 1). The column structure is determined by a set of 
equations including local mass conservation law and hydro¬ 
static equilibrium. The mass conservation law can be ex¬ 
pressed as follows: 


M 

25d 


= pv, 


(15) 


where p and v are local mass density and velocity cor¬ 
respondingly. The exact local values of p and v cannot 
be calculated without solving the full set of radiation- 
hydrodynamic equations. The velocity profiles in a sinking 
region were calculated for particular cases (Wang & Frank 
1981) and can be approximated by power law v « with 
the exponent £ ~ 1 for relatively low accretion luminos¬ 
ity ~ 10 3 ' ergs -1 and with higher £ € [1; 5] depending on 
the accretion channel thickness, opacity in a sinking region 
and magnetic field structure. For the case of high accretion 
column (H ~ 3.5 R) and a = ot the velocity profile was cal¬ 
culated to be u « h 5 (Basko & Sunyaev 1976). Thus, we as¬ 
sume that the velocity profile in the accretion column below 
the shock zone is given by a power law and discuss influence 
of the parameter £ on the solutions (see Section 3.2). The 
velocity at the NS surface is equal to zero, while the plasma 
velocity at the top of the column just below the shock is 
v(H) = vg(H)/7 (see Lyubarskii & Syunyaev 1982; Becker 
1998), where vs{H) = \J2GM/(R -I- H) is a free-fall veloc¬ 
ity at height H above the NS surface. 

According to equation (1), the local radiation pressure 
Trad can be presented as follows: 


Trad ( t X 1 h) 


f GM 2 F EAA (H X ) 

/'isw d ’b— 


(16) 


where we used standard boundary condition for the radia¬ 
tion density e ra d = 2 F/c and the approximation of thermo¬ 
dynamic equilibrium (12). The connection between radiation 
pressure P lad at the centre of the accretion channel and the 
emergent radiation flux Tj_,esc is determined approximately 
by the radiative transfer equation (4). However, for the accu¬ 
rate solution one has to account that the opacity k± depends 


on the local column conditions such as the temperature and 
the B-field strength (see details in Section 2.3). Equation (4) 
can be transformed to 


Trad (•£, fi) 


2 p Tj_,esc(T) 
dh C 


d-h/2 

/ 


1 2 F_ L esc 

K± zdz-h- --— 

3 c 




(17) 


where the similar surface boundary condition like in equa¬ 
tion (16) was used. Because k± depends on the local tem¬ 
perature which is evaluated from P lad (x,h) (see eq. (12)), 
equation (17) could be solved iteratively, if one knows the 
radiation pressure P ra d at the centre of the accretion chan¬ 
nel. Then the flux emergent from the wall of the column 
Tj_,esc can be found for the known column thickness dh and 
plasma mass density p at every height. The thickness and the 
mass density are found iteratively as well (see the iterative 
scheme below). 

The described model could be incorrect for the case of 
extremely heigh accretion columns (H R) due to geomet¬ 
rical reason, namely due to curvature of dipole magnetic field 
lines, when the magnetic field and gravity are not aligned 
and condition given by eq. (1) is not valid. It means that 
we cannot find any solution at some high luminosities when 
H > R. Fortunately, the approximate formulae (8) and (9) 
show that the luminosity saturates at H as R and we still 
can obtain a sufficiently good estimation of the maximum 
X-ray pulsar luminosity for given NS mass and radius and 
the magnetic field strength. 


2.3 Scattering cross-section 

In order to compute the accretion column luminosity or the 
accretion column height for a given luminosity, we need to 
find the mean opacity across k± and along km the magnetic 
field. We assume that the plasma in the column is fully ion¬ 
ized and ignore bremsstrahlung opacity. Thus, the plasma 
opacity is determined by Compton scattering only. Comp¬ 
ton scattering cross-section in case of high magnetic field de¬ 
pends strongly on the photon energy, polarization state and 
photon momentum direction (Daugherty & Harding 1986; 
Mushtukov, Nagirner & Poutanen 2015). In the case of op¬ 
tically thick plasma, in diffusion approximation, we only re¬ 
quire the Rosseland mean opacities in the two directions: 


_ n B 
KX,|| = 0 -_l.ii — 


o--i-.ll 

p e mn ’ 


(18) 


where cTl,|| are the Rosseland mean Compton scattering 
cross-sections across and along magnetic field (see below), 
p e = 1.17 is the mean number of nucleons per electron for 
fully ionized solar mix plasma. For the angular-dependent 
case, the cross-section along and across the 13-field are 
(Paczynski 1992; van Putten et al. 2013; Mushtukov et al. 
2015) 


OO 

/ 


dB E (T) 

dT 


dE J dp 3 p 2 


aj{E,p,T) 


f 

Jo 


d B e (T) 
dT 


(19) 


dE 
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and 


oo 7r 


dE dip d/in-Mn 


7T n (jj(E,n,T) 


cr_L, 


/■ 

Jo 


d Be 
d T 


; ( 20 ) 


d E 


respectively. Here T is the local electron temperature, Be(T ) 
is the Planck function, E is a photon energy, /i is a cosine 
of the angle between photon momentum and the B-held di¬ 
rection and is a cosine between the momentum and per¬ 
pendicular to the held direction (/r = y/l — cos ip). Index 
j correspond to a given photon polarization state (j = 1 for 
X- and j = 2 for O-mode). In case of mixed polarization the 
effective cross-section can be estimated as follows 


l/a = f/ai + {1- f)/a 2 , ( 21 ) 

where / is a fraction of radiation in the X-mode 
(Mushtukov et al. 2015). We assume / = 1 for further 
calculations. It is a rather good approximation because 
most of the radiation is expected to be in X-mode due 
to the effective swapping of photons from O- to X-mode 
(Miller 1995; Mushtukov et al. 2012). In case of sufficiently 
strong magnetic held the typical photon energy E can be 
much lower than the cyclotron one E cyc \ ~ 11 . 6 Bi 2 keV. 
In case of E < E cyc i the scattering cross-sections can 
be approximated as cr±(E) ~ ctt(E / E cyc \) 2 and <J\\{E) ~ 
ot[(1 — f 2 ) + lT(E/E cyc \) 2 ] for X- and O-mode respectively 
(Canuto et al. 1971). As a result the effective cross-section 
for extremely strong S-held tends to be signihcantly lower 
that the Thomson value. 

The computed values of ct(_l , 11 ) ( 1 , 2 ) were tabulated for 
40 values of the magnetic held strength distributed loga¬ 
rithmically in the range 10 11 — 10 15 G. The opacities were 
calculated for 300 values of the temperature in the range 
10 6 — 10 9 K for every magnetic held strength B. The vari¬ 
able steps were specially adopted for more detailed descrip¬ 
tion in the temperature region near corresponding electron 
cyclotron energy. 


2.4 Accretion column base geometry 

It is assumed that the accretion how is interrupted at 
the magnetospheric radius due to the sufficiently strong 
magnetic held pressure. The inner radius of the accretion 
disc depends on the magnetic held strength and its struc¬ 
ture, mass accretion rate and the feeding mechanism of the 
pulsar. It can be estimated with the following expression 
(Ghosh & Lamb 1978, 1979, 1992) 

R m = 7 x 10 7 Am 1 / 7 i?g 0 / 7 Bi 2 / 7 L 3 9 2/7 cm, ( 22 ) 

where A is a constant depending on the accretion how geom¬ 
etry (A = 1 for the case of spherical accretion and A < 1 for 
the case of accretion through the disc, with A = 0.5 being 
a commonly used value), m = Mj Mg and we use notations 
in cgs units for other variables Q = Q x 10 x . 

Accretion channel geometry is dehned by the interac¬ 
tion between the accretion disc and NS magnetosphere. The 
thickness of the accretion channel depends on the penetra¬ 
tion depth 5 of the accretion disc into the NS magnetosphere. 
It is expected to be of the order of the disc scale-height at 
the inner edge Hd- The accretion disc structure at the inner 


edge depends on the B-held strength and the luminosity. 
We evaluate Hd using the Shakura & Sunyaev (1973) model 
with a slightly different vertical averaging and with the ac¬ 
curate Kramer opacity (Suleimanov, Lipunova, & Shakura 
2007). The disc can be divided into three zones according 
to the dominating pressure and opacity sources. The gas 
pressure and the Kramer opacity dominate in the outer re¬ 
gions (C-zone), the gas pressure and the electron scattering 
dominate in the intermediate zone (B-zone) and the radia¬ 
tion pressure may dominate in the inner zone (A-zone). The 
boundaries between A- and B- and between B- and C- zones 
are expected to be at tab ~ 2.2 x 10 8 L^ 21 ' Rg 6 ^ 21 cm 

and rec ~9x 10 8 Z/ 2 9 3 r?W 1 / 3 .Rg ' /3 cm, respectively. The disc 
structure at the inner edge is dehned then by the magneto- 
spheric radius f? m . 

A relative disc scale-height at radius r for the C-zone, 
where accretion disc is expected to be interrupted in case of 
relatively low mass accretion rate, is 

« 0.056 a " 1/10 L% 20 m- 21/4 ° J? 3/2 ° r 4/s . (23) 

c 

Assuming that matter leaves the disc at magnetosphere ra¬ 
dius and follows the dipole magnetic held lines down to the 
polar caps of the NS we get the ratio lo/do and the area Sb- 
For the case of a truncation in the C-zone: 

(If) «2tt^.~95A ~ 1 / 8 m 71 / 140 B 12 1 / 14 L 3 9 4/35 , (24) 

\ do J c H a 

and corresponding area of the accretion channel base 
st C) « 6.6 x 10 9 A _ 7 / 8 L 2 9 5 B(” 2 1 / 2 m_ 13 / 20 ^ 6 9/10 c m 2 - (25) 
A relative disc scale-height in a relatively narrow B-zone is 

* 0.064 a - 1/10 Lg 9 5 m- 11/2 ° J ? 4/5 r 4/2 °, (26) 

B 

and the ratio lo/do can be estimated as follows 

(J^j «80A ■ 1 / 20 m 19 / 3 S B 1 " 2 1 / 35 f 39 13/70 . (27) 

The accretion channel base area in this case is 

4 B) ~ 7 x 10 9 A- 19 / 20 L 4 ®/ 70 S- 41 / 70 m - 24/35 ^ 29/70 cm 2 . 

(28) 

In the case of sufficiently high luminosity, L 39 )> 

0.34 A 21/,Z2 the structure of the accretion 

disc zone A is dehned mostly by the radiation pres¬ 
sure and the disc scale-height is independent of radius: 
(Shakura & Sunyaev 1973; Suleimanov et al. 2007): 

u ~ 3 M 7 T 39 -R 6 , 00 n 

H d ,A «kt- - ^ 10 - cm (29) 

87 r c m 

and therefore 

( J) a ~ 2tt^ ~ 44 A m 8/7 R 3 6 /7 Bti 7 L^ 9/T . (30) 

The corresponding area of the accretion channel base 

st A) = 1.3 x 10 10 A~ 2 Lll /7 B^ 2 8/7 m- 9/7 Rl /7 cm 2 . (31) 

If the magnetic dipole is inclined with respect to the 
orbital plane, the expression for the hot spot area and the 
ratio (lo/do) would be different and the spot would have a 
shape of an open ring. We use an additional parameter Iq/Io, 
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which shows what part of the full ring length lo is filled with 
accreting gas. In numerical calculations we take Iq/Io > 0.5. 
It is also natural to assume that the ratio Iq/Io is higher for 
higher mass accretion rates. 


2.5 Computation scheme 

A general computation scheme consists of a few steps. The 
main task is to find the central accretion column height H, 
which corresponds to a given X-ray pulsar luminosity L (i.e. 
to a given mass accretion rate M). The parameters of the 
model are NS mass M and radius R, surface magnetic field 
strength B, parameter £, which determines the velocity pro¬ 
file in the sinking region, the polarization composition of 
radiation given by parameter / (see Section 2.3) and the 
ratio Iq/Io, which gives what part of the full ring is filled 
with the accreting matter (see Section 2.4). The accretion 
column base geometry (lo, do) is defined by L and B (see 
Section 2.4). The accretion column height and the structure 
are calculated by double iteration scheme. The external it¬ 
eration finds the accretion column height for a given set 
of parameters by the dichotomy method, while the internal 
iteration calculates the luminosity for current accretion col¬ 
umn height and given base geometry. We make a first guess 
of the accretion column height using approximation (8) and 
assuming k± = kt. The steps of the iterative scheme are: 

(i) The temperature inside the sinking region T(x,h) is 
calculated analytically using equations (3), (5), (6) and (12) 
as a first guess and recalculated iteratively further (see be¬ 
low). 

(ii) For the current temperature distribution we compute 
the opacity k,\\,±( x,h) (see Section2.3). 

(iii) For the current accretion column height H we com¬ 
pute the central radiation pressure (x = 0) using formula 
(16) as a function of height h inside the sinking region, i.e. 
we get P ra d(0, h) for h G [0;id]. Formally, formulae (15)- 
(16) lead to the infinite values of the mass density and the 
radiation pressure at h —► 0. It means that formula (16) is 
inaccurate at the very bottom of the accretion column. The 
total pressure becomes higher than the magnetic pressure 
(Trad + Tg) > B 2 /4 -7T and plasma spreads over the NS sur¬ 
face. The height ho, where it happens, is taken as the base 
of the accretion column. 2 

(iv) Then we find a dependence of the thickness of a 
sinking region dh on height h. For that, we first compute 
FT,esc (h = ho) using formula (17) from the already known 
Trad(0, ho), current opacity k±(x, h) and known thickness at 
the base dh 0 = do- Then using the same formula (17) and 
already known F± tesc (h = ho) we compute the radiation 
pressure at the base h = ho from the middle of the column 
to the edge, i.e. we get P ra d(®, h = ho) for x G [0; do/2]. The 
evaluated radiation pressure at the accretion column base 
Trad (x, h = ho) is used to get the column height H x at dis¬ 
tance x from the center of the channel (see equation (16)). 
Then we get P v a d(®, h) using again equation (16). Reversing 


2 The height ho depends on the mass accretion rate, B-field 
strength and velocity profile in the sinking region and takes on 
a value ~ ( 0.1 — 100 ) cm, which is -C do and does not affect the 
luminosity calculated for given accretion column height. 


the function H x we finally get the dependence of a sinking 
region width dh on the height h. 

(v) As soon as the central radiation pressure P T ad(0,h) 
and the thickness of the sinking region dh are known at 
any height, we use equation (17) and get the emergent flux 
F±,esc(h) as a function of h. 

(vi) We get the resulting X-ray pulsar luminosity for a 
given accretion column height H and the temperature struc¬ 
ture 

H 3/2 

L' ~4Z 0 J F±Mh)dh. (32) 

o 

(vii) Using computed P IS ,d(x,h) we get new temper¬ 
ature distribution in the accretion column T(x,h) = 
(3F ra d(i,/i)/tt)^ 4 . Then we return to step (ii), recalculate 
the luminosity L' for a given column height H and new tem¬ 
perature T(x, h) and continue this procedure until we get the 
accretion luminosity L'(H) with the relative accuracy of 1 
per cent. 

(viii) If the luminosity L'(H) does not coincide with the 
given accretion luminosity L we get the next guess of the 
accretion column height and return to step (i). We continue 
the iteration procedure until the relative accuracy of 1 per 
cent is reached. 

Finally, we get the height and the structure of the accretion 
column for a given mass accretion rate. 


3 NUMERICAL RESULTS 

3.1 Accretion column properties 

There are four parameters (in addition to the NS mass M 
and radius R) defining the accretion column height and 
structure. These are the surface polar NS magnetic field 
strength B, the mass accretion rate M or the accretion lu¬ 
minosity L = GMM/ R, and the geometrical sizes of the 
accretion column base, do and lo- Here we investigate how 
these parameters affect the properties of the accretion col¬ 
umn. In order to separate the effects caused by different 
reasons, we will assume a fixed accretion column base ge¬ 
ometry with the following parameters lo = 7.7 x 10 5 cm 
do = 1.3 x 10 4 cm in this section only. The self-consistent 
results, where the dependence of accretion channel geome¬ 
try on the mass accretion rate is taken into account, will be 
discussed in Section 3.2. 

Two main parameters that determine the global proper¬ 
ties of the accretion column are the local mass accretion rate 
M/2Sd and the effective scattering cross-section across the 
accretion channel a±. The vertical temperature structure in 
the centre of sinking region of the accretion column for a 
given column height depends on the local mass accretion 
rate only. The temperature structure across the accretion 
channel and, therefore, the emergent flux at given height is 
determined by the local effective scattering cross-sections 
ctj_(i 2 ) (see Section 2.3) and the column density pdz 
across the channel. The stronger the field, the smaller the 
scattering cross-section and the smaller the optical thick¬ 
ness across the column. As a result, the higher the R-field, 
the higher the emergent flux F± for the fixed temperature in 
the column centre and the column density. Global properties 
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Figure 2. Theoretical dependences of (a) the internal tempera¬ 
ture, (b) the effective surface temperature, and (c) the effective 
scattering cross-section on the h/R ratio for L = 10 39 ergs 1 
and fixed accretion column base geometry: Iq = 7.7 X 10 s cm 
do = 1.3 X 10 4 cm. Red solid, green dotted, and blue dashed lines 
correspond to B = 3 X 10 12 , 10 13 and 3 X 10 13 G, respectively. 
Parameters: M = 1.4Mq, R = 10 6 cm, £ = 1, / = 1. 


of a self-consistent accretion column model, which have to 
be constructed in an iterative way, are determined by these 
facts. 

Let us focus on the dependence of accretion column 
properties on the magnetic field strength and consider ac¬ 
cretion column models corresponding to three different sur¬ 
face B-field strength and fixed accretion luminosity L = 
10 39 ergs -1 (see Fig. 2). The higher field strength leads to 
lower effective scattering cross-section (Fig. 2 c). Therefore, 
the emergent flux at a given height (F± = assT/lg) is higher 
(Fig. 2 b). Smaller column at high B-field provides the same 


luminosity as higher column at relatively low field strength. 
On the other hand, the increased leaking of thermal energy 
makes the central column temperature lower in case of high 
B-field (Fig. 2 a) and supporting of the accretion column is 
less effective. Then the accretion column is also smaller in 
case of higher field strength for fixed mass accretion rate. 
We note, that the significant dependence of accretion col¬ 
umn properties on the magnetic field strength takes place for 
sufficiently high magnetic field strengths (B > 10 12 G) only, 
when most energy is radiated below the electron cyclotron 
energy E cyc = 11.6 B \2 keV and the effective electron scat¬ 
tering cross-section drops significantly below or (Fig. 2 c). 
As B < 10 12 G, the properties of sufficiently bright accre¬ 
tion columns become almost independent on the magnetic 
field strength. 

The effective temperature of the accretion column is not 
a monotonic function of height and reaches its maximum 
value above the NS surface (Fig. 2b). On the contrary, the 
internal temperature decreases monotonically with height 
(Fig. 2a), according to the basic principles described above. 
The value of internal temperature on the top of the column 
is defined by the boundary condition in equation (16). The 
internal central temperature is higher than the effective tem¬ 
perature at a given height and this difference is defined by 
the optical thickness across the accretion column. The op¬ 
tical depth across the column decreases with the height be¬ 
cause of decreasing of geometrical thickness (see Section 2.1) 
and overcompensates the drop of the internal temperature. 
This causes the increase of the effective temperature with 
height. 

We note that the difference between the effective and 
the internal temperatures is larger for a relatively low B- 
field strength at a fixed luminosity, because of higher elec¬ 
tron scattering opacity. These temperatures become close to 
each other as soon as the optical thickness across the column 
becomes close to unity. It happens when the accretion lumi¬ 
nosity is close to the critical luminosity value L* for a given 
magnetic field strength (Mushtukov et al. 2015), which de¬ 
fines the minimum luminosity when the accretion column 
exist. 

It was shown above that the column reaches its maxi¬ 
mum luminosity when the column height H ~ R. Therefore, 
a stronger magnetic field allows X-ray pulsars to reach higher 
luminosities, i.e. the accretion column of a fixed height will 
be more luminous when the magnetic field is stronger. We 
illustrate this statement by a few accretion columns with 
the same accretion column base geometry and different lu¬ 
minosities, for two fixed magnetic field strengths: 10 13 and 
10 14 G (see Fig. 3). It is clear that the accretion columns of 
similar heights but different field strength produce luminosi¬ 
ties that may differ by an order of magnitude. 

The dependences of column properties on the geomet¬ 
rical parameters are much simpler and we do not illustrate 
it in details. Decrease of the column base thickness do and 
the increase of column base length lo lead to lower accre¬ 
tion column heights at a given mass accretion rate and the 
B-field strength. In fact, the geometrical parameters are not 
independent, they depend on the magnetospheric radius Rm 
and the accretion disc scale-height at its inner edge. Both of 
them depend on the mass accretion rate and the magnetic 
field strength on the NS surface. We consider this problem 
in detail in the next section. 
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Figure 3. The internal temperature and the effective scattering cross-section as a function of height. Different curves correspond to 
various accretion column height: 0.25 R, 0.5 R, R and different luminosities, which depend on the field strength. The accretion column 
base geometry is fixed: lo = 7.7 X fO 5 cm, do = 1.3 X 10 4 cm. Red solid lines correspond to the case of surface magnetic field strength 
B = 10 13 G, while blue dashed lines correspond to B = 10 14 G. The corresponding luminosities (in units of 10 39 ergs ') are marked 
next to each curve. The higher the accretion column, the higher the internal temperature, the higher effective scattering cross-section 
for a given height above the NS surface. Parameters: M = 1.4Mq, R = 10® cm, £ = 1, / m 1. 



H/R 

Figure 4. The accretion column luminosity as a function of its 
height. Red solid lines correspond to various surface magnetic field 
strength: B = 10 11 ,3X10 11 ,10 12 ,3X10 12 ,10 13 ,3X10 13 ,10 14 ,3X 
10 14 and 10 15 G. The red circles at the right end of the curves 
indicate the maximum accretion luminosity values and the corre¬ 
sponding accretion column height. The red squares indicate the 
minimum critical luminosity L* when the accretion column ex¬ 
ists. The saturation of the luminosity is caused by two different 
factors: a switch of the accretion disc behaviour on its inner ra¬ 
dius for B < 10 13 G and by increasing of the effective scattering 
cross-section for B > 10 13 G, which makes the column growth 
ineffective at some point. The black dotted line shows the approx¬ 
imate dependence (9) of the luminosity on the accretion column 
height for k± = kj. The parameters used in the calculations are 
M = 1.4M 0 , R = 10 6 cm, A = 0.5, l/,/l 0 = 1, £ = 1, / = 1. 

3.2 Maximum luminosity 

Using the accretion column model described above we inves¬ 
tigate the maximum accretion luminosity of magnetized NS 
as a function of the surface magnetic field strength. Higher 
accretion luminosity corresponds to a taller accretion col¬ 
umn. 


The velocity profile inside the sinking region, which is 
defined by the parameter £, does not affect much the lumi¬ 
nosity dependence on the accretion column height: its vari¬ 
ations within the interval 1 < £ < 5 lead to variations of 
accretion luminosity by a factor <( 2 for a given column 
height. The higher the £, the lower the calculated lumi¬ 
nosity value. If one takes scaled velocity structure obtained 
by Basko & Sunyaev (1976) for the accretion column with 
constant scattering cross-section a = ot, the luminosity is 
smaller by a factor < 2 than that given by £ = 1. Thus, we 
take £ = 1 in our calculations. 

Theoretical dependences of X-ray pulsar luminosities on 
the accretion column height for various 13-field strength are 
shown in Fig. 4. All dependences start at relatively small 
accretion column heights corresponding to the critical lu¬ 
minosity L*, which depends on the 13-field strength as well 
(Mushtukov et al. 2015). The lowest height is expected to 
be comparable with the thickness do of the accretion chan¬ 
nel near the NS surface. Below the critical luminosity, the 
energy is emitted by a hotspots at the NS surface. The lu¬ 
minosity dependences on the height of the accretion column 
are close to the approximate analytical dependence L(H/R) 
given by equation (9) for the case of relatively low magnetic 
field strength (B < 3 x 10 12 G). The deviations of the com¬ 
puted dependences from the analytical one can be explained 
by changes of accretion channel geometry, which is given 
by the ratio lo/do assumed to be constant in the approxi¬ 
mate relation (9). The dependence of lo/do on the accretion 
rate varies according to the region where the accretion disc 
is interrupted by the NS magnetic field (see Section 2.4). 
If the accretion rate is relatively low and the disc is in¬ 
terrupted in the C-zone, the accretion column luminosity 
L oc f 33 ^ 39 (H/R). However, if M becomes sufficiently high 
and the disc is disrupted in the A-zone, the luminosity grows 
with the column height much slower as L oc f 7 ' 16 (H/ R). 
This transition from the C- and B-zones to the A-zone is 
clearly seen in Fig. 4. The transition luminosity depends on 
the magnetic field strength and is given by L% ~ 0 .2 B^ 11 . 

The luminosity dependences on the accretion column 


MNRAS 000, 1-11 (2015) 













The maximum luminosity for magnetized neutron stars 9 


height L(H) for B > 10 13 G are shifted to higher values 
because of a strong reduction of the effective electron scat¬ 
tering cross-section. The averaged scattering cross-section 
falls quickly at high B and, therefore, the factor kt/k± 
grows substantially. We see, however, that the luminosity 
does not grow as much with the column height as it hap¬ 
pens at lower B. This is related to the fact that with in¬ 
creasing M the optical thickness across the column grows 
leading to a higher inner temperature and correspondingly 
higher scattering cross-section (because the photon peak be¬ 
comes closer to the resonance) which leads to a reduction of 
the escaping flux. Thus at some point the further growth of 
the column is not effective any more and luminosities higher 
than some maximum luminosity (for given B-field strength) 
cannot be supported by accretion column of any height. 

The dependence of the maximum luminosity on the sur¬ 
face magnetic field is shown in Fig. 5 by a thick solid line. 
We limited the maximal accretion luminosity for low S-field 
strength by hand at H/R = 1, because of a much slower lu¬ 
minosity grow at H/R > 1 and a questionable validity of our 
model for such high columns. However, in any case we see 
that the maximum column luminosity L^g ~ 0 . 35 . 8^2 4 can 
easily reach value in excess of 10 4 ° ergs -1 for magnetar-like 
magnetic fields. The parameters of the model such as the 
exact velocity profile in the sinking region of the accretion 
column and the fraction of photon in each polarization state 
do not affect much the final results. The relation between the 
luminosity and the B-field can be reversed implying that for 
a given luminosity (L39 > 2 ) there is a lower limit on the 

magnetic field: 812 )> . 

4 DISCUSSION 

4.1 An admissible accretion luminosity range 

We have shown in the previous section that a tall accretion 
column can radiate in excess of 10 4 ° ergs -1 for the magnetar 
like magnetic fields B 10 14 G. However, we have not ad¬ 
dressed the question if a necessary high accretion rate is con¬ 
sistent with the presence of the column. The problem is that 
at high accretion rate the accretion disc may become locally 
super-Eddington outside of the magnetosphere producing 
strong outflows (Shakura & Sunyaev 1973; Poutanen et al. 
2007) and reducing the NS luminosity. Additional problem 
arises because in this case the disc becomes very thick at 8 m 
leading to a nearly spherical accretion, which is limited by 
the Eddington luminosity. The correction of the opacity at 
high B will not play any role here just because the magnetic 
field drops rapidly with the distance and at some radius the 
cyclotron resonance will appear in the X-ray range where 
most of the photons are escaping. 

The condition for the absence of the outflows outside 
the magnetosphere can be written as 

R sp < Rm, (33) 

where R sp is the spherization radius where the radiation 
force is still balanced by gravitation. Its simplest estimate 
is based on equating the scale-height of the radiation- 
dominated disc to the radial distance (Shakura & Sunyaev 
1973): 

R sp = — — M = 7.5 x 10 6 L39Re cm. (34) 
8n c m 



Figure 5. The maximum accretion luminosity L** as a func¬ 
tion of the magnetic field strength is given by thick solid line. 
In the interval 10 13 < B < 10 15 G it can be approximated as 
LJg ~ 0.35B^2 4 - The lower limit on the X-ray luminosity (37) 
shown by dashed lines is related to the inhibition of accretion 
by the propeller effect. The upper limit (35) shown by the dot¬ 
ted line corresponds to the condition that the spherization radius 
is smaller the magnetospheric radius. The critical luminosity L * 
when the accretion column starts to grow is shown by the dashed- 
dotted line (Mushtukov et al. 2015). The most likely position of 
the ULX in galaxy M82 (Bachetti et al. 2014) is given by the 
orange circle, which corresponds to B ~ 8 x 10 13 G, H/R ~ 1, 
lo = 4 x 10 5 cm and do — 1.7 x 10 4 cm. Parameters: M = 1.4Mq, 
R = 10 6 cm, A = 0.5, Iq/Io = 1, £ = 1, / = 1. 

It may be smaller by a factor of a few if one accounts for the 
energy advection (Poutanen et al. 2007). Using the condi¬ 
tion (33) we get an approximate upper limit for the accretion 
luminosity: 

L™ ~ 5.7 A 7/9 B 4 ^ 9 m 8/9 i? 1/3 (35) 

shown by a dotted line in Fig. 5. We see that the theoret¬ 
ical maximum column luminosity is always below the limit 
coming from this spherization constraint and therefore our 
model is self-consistent. 

In addition to the upper limit on the luminosity, there 
should exist a lower limit too. At low accretion rates, the 
spin frequency of the NS may exceed the Keplerian fre¬ 
quency at the magnetospheric radius. Then the accretion 
of matter may be inhibited by the centrifugal force exerted 
by the rotating magnetosphere and the matter then may 
be expelled from the system due to the “propeller effect” 
(Illarionov & Sunyaev 1975). This dramatically reduces the 
accretion luminosity. The propeller starts to operate when 
the magnetospheric radius R m becomes larger than the coro¬ 
tation radius 

(GMP 2 \ ! s 1/3 d 2/3 

Rc = I — j — 1 ~ 1.5 x 10 m ' P ' cm, (36) 

where P is a NS spin period in seconds. This happens at 
luminosity 

i“ 9 in ~ 7 x 10 -2 A 7/2 m~ 2/3 R 5 6 Bf 2 P~ 7/3 , (37) 

which defines the likely lower limit on X-ray luminosity of a 
binary system. 
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During the accretion process a NS in a binary sys¬ 
tem is spun up to the period determined by the condi¬ 
tion R m ~ R c , when the accretion luminosity is given by 
equation (37). For the Eddington accretion rate relaxation 
happens very rapidly, within r re 1 ~ 300 1 15 m} ^ 3 3 yr, 

where 745 = 7/10 45 gem 2 is the NS moment of inertia and 
7^30 = /r/10 30 Gem 3 is the NS dipole magnetic moment 
(Lipunov 1982). Of course, the luminosity may exceed the 
minimum value (37) over short time intervals when the ac¬ 
cretion rate increases. 

We conclude that the X-ray pulsar luminosity is likely 
to be bound from below by L mm and from above by L** « 
0 . 3573 ^ 2 - ^ is clear from Fig. 5 that the stronger the mag¬ 
netic field, the narrower the possible luminosity range for a 
given NS spin period. A likely maximum accretion luminos¬ 
ity given by the interception of the two constraints is then 

L 39 ss 1.6 m 2 / 5 P 5 / 4 i?g 3 A -21 / 10 , (38) 

corresponding to the magnetic field 

B 12 m 7.5 m 8 / 15 P 28 / 15 i?g 4 A -14 / 5 . (39) 

Thus the luminosities of X-ray pulsars in excess of 
10 4 ° ergs -1 would definitely require an extremely strong 73- 
field (73 > 10 14 G) and a large spin period (P > Is). 

4.2 X-ray pulsars and ULXs 

All known X-ray pulsars can roughly be divided into two big 
groups: persistent and transient sources. Among the bright¬ 
est persistent sources we mention here Cen X-3, SMC X-l 
and LMC X-4. The only source in this group with the known 
magnetic field is Cen X-3, which shows the cyclotron line 
around 28 keV (Nagase et al. 1992) that corresponds to the 
corrected for the redshift magnetic field 73 ~ 3.2 x 10 12 G. 
The brightest X-ray pulsars in the Magellanic Clouds do not 
show any cyclotron features in their spectra and hence the 
magnetic fields there are unknown. The maximal observed 
luminosities in this group range from ~ 5 x 10 3 ' ergs -1 for 
Cen X-3 to - 10 39 ergs -1 for SMC X-l and LMC X-4 during 
outbursts. 

Another group of transient sources mainly consists of 
NS in binary systems with the Be companions. They are 
known to exhibit rare non-periodic type II outbursts with 
the peak X-ray luminosity reaching the Eddington value (for 
a review see Reig 2011). The brightest sources in this group 
are 4U 0115+63 and V 0332+53 with the peak luminosi¬ 
ties as high as a few xl 0 38 ergs -1 and a clear evidence of 
the high accretion column (see e.g., Tsygankov et al. 2006, 
2007; Tsygankov, Lutovinov & Serber 2010; Poutanen et al. 
2013; Postnov et al. 2015). The magnetic fields are known 
from the cyclotron lines: ~ 1.3 x 10 12 G for 4U 0115+63 
(White, Swank & Holt 1983) and ~ 3.2 x 10 12 G for V 
0332+53 (Makishima et al. 1990). 

Apart of these two groups of X-ray pulsars we can men¬ 
tion a unique transient source GRO J1744-28 showing both 
X-ray pulsations and unusual short bursts. The broadband 
observations of this source revealed a cyclotron absorption 
line in the energy spectrum at 73 ~ 4.5 keV (D’Ai et al. 
2015; Doroshenko et al. 2015), which corresponds to 73 ~ 
5 x 10 11 G. The maximum luminosity during the outbursts 
was ~ (3 + 4) x 10 38 ergs -1 . The observed luminosities of 
all these pulsars are in good agreement with our limits. 


It was shown above that the accreting NS may have 
accretion columns radiating in excess of 10 4 O ergs -1 for 
magnetar-like magnetic fields Bo !> 10 14 G (see Fig. 5). 
In other words, some ULXs can actually be extreme X-ray 
pulsars as indeed was recently observed for X-2 source in 
galaxy M82 (Bachetti et al. 2014). Its high pulsating lumi¬ 
nosity ~ 4.9 x 10 39 ergs -1 implies according to our model 
the magnetic field 73 > 10 14 G. The observed pulsation pe¬ 
riod of 1.37 s makes the magnetic field of ~ 10 14 G the only 
possible option for this particular ULX, because the higher 
field strength leads to the propeller effect (see Section 4.1). 
Because the ULX is very likely to be close to the equilibrium 
state (Bachetti et al. 2014; Lipunov 1982), an estimation of 
its magnetic field based on the analysis of the period P and 
the period derivative P is not valid since one need to know 
exact physics of disc and magnetosphere interaction and 
the estimations of magnetic field strength become strongly 
model dependent near the equilibrium (Parfrey et al. 2015; 
Dall’Osso et al. 2015; Ek§i et al. 2015). Transition between 
accretion and propeller regime has been detected recently 
in M82 X-2 (Tsygankov et al. 2015) confirming our estima¬ 
tion of the magnetic field strength of ~ 10 14 G. A high ac¬ 
cretion rate needed to power high ULX luminosities (L ~ 
10 39 — 10 41 ergs -1 ) from pulsars imply a large thickness of 
the accretion disc on its inner edge. This results in repro¬ 
cessing of a significant part of pulsar radiation in the disc, 
which might lead to spectral deviations comparing to the 
intrinsic spectrum of the accretion column. 


5 SUMMARY 

We have constructed a simplified model of the accre¬ 
tion column on a magnetized NS. We based the model 
on previous suggestions by Basko & Sunyaev (1976) and 
Lyubarskii & Syunyaev (1988) that the column is vertically 
supported by the Eddington flux. We use diffusion approxi¬ 
mation to evaluate the radiation flux through the sides of the 
column. For the first time, we use the exact Compton scat¬ 
tering cross-sections in strong magnetic field and find the 
column temperature structure and the cross-sections self- 
consistently via an iterative procedure. As a result we have 
obtained the accretion column luminosity as a function of 
the magnetic field strength and the column height. We also 
evaluate for a given surface magnetic field the maximum pos¬ 
sible column luminosity that is reached when the column 
height becomes comparable to the NS radius. We showed 
that the column may radiative in excess of 10 4 O ergs -1 in 
the case of the magnetar-line magnetic field > 10 14 G. 

We also evaluate the possible range of accretion lumi¬ 
nosity in X-ray pulsars. It is limited from above by the con¬ 
dition that the accretion disc does not spherize outside the 
magnetosphere. This condition, however, gives the limiting 
luminosity above our upper limit thus confirming the con¬ 
sistency of our calculations. The pulsar luminosity may also 
be limited from below by the propeller effect. 

The range of possible pulsar luminosities becomes nar¬ 
rower for stronger magnetic fields. The surface magnetic field 
strength can be limited from below if one knows the lumi¬ 
nosity of the object: 73 i 2 )> 4 Egg 3 . Thus the X-ray pulsar lu¬ 
minosity exceeding a few 10 39 ergs -1 requires magnetic field 
> 10 13 G, while the ULX luminosities of 10 4 °ergs -1 would 
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be consistent only with the magnetar-like fields > 10 14 G. 
On the other hand, the magnetic field cannot be too large, 
because of the propeller effect. We showed that the estimated 
X-ray luminosity of ULX X-2 in galaxy M82 (Bachetti et al. 
2014) of a: 10 40 ergs _1 can be consistent with our model 
only for B « 10 14 G. 

Because of the set of constraints it is likely that the 
luminosity L ~ (1 4- 3) x 10 40 ergs _1 is a good guess 
for the NS maximum possible accretion luminosity. Inter¬ 
estingly, this value coincides with the cut-off observed in 
the X-ray luminosity function of high-mass X-ray binaries, 
which does not show any features at lower luminosities 
(Mineo, Gilfanov & Sunyaev 2012). Therefore, it is natural 
to suggest that accreting NSs are not rare among ULXs and 
their fraction is probably close to their fraction among the 
high-mass X-ray binaries. 
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